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INTRODUCTION 


The  shelf  life  of  ammunition  is  affected  primarily  by  temperature,  relative  humidity,  and 
solar  radiation.  Temperature  control  is  essential  for  maintaining  ammunition  serviceability. 
Characterizing  storage  facilities  to  optimize  magazine  usage  and  design  can  be  accomplished 
by  maximizing  insulation  and  proper  airflow,  thereby  improving  temperature  control.  Equation 
based  modeling  will  allow  storage  facilities  to  be  quantitatively  characterized  and  thereby  opti¬ 
mized  for  maintaining  high  quality  ammunition  and  extending  shelf  life. 

Scope 

This  study  will: 

•  Quantify  the  adverse  environmental  effects  on  ammunition  in  storage 

•  Derive  engineering  equations  to  determine  the  temperature  conditions  that 
ammunition  is  exposed  to  within  storage  magazines 

•  Develop  life  prediction  models  to  forecast  ammunition  shelf  life 

Objectives: 

•  Develop  environmental  equations  that  delineate  inside  storage  temperature  at 
various  geographic  locations  as  a  function  of  diurnal  ambient  temperature, 
wind  velocity,  and  solar  radiation 

•  Demonstrate,  by  example,  the  integration  of  the  developed  environmental 
equations  with  existing  life  prediction  models  to  forecast  shelf  life 

Approach 

The  goal  is  to  develop  models  that  predict  ammunition  shelf  life  based  on  the  storage 
facility,  location,  duration  of  storage,  and  previous  stockpile  exposure.  To  accomplish  the  goal, 
life  models  were  developed  based  on  historical  data.  Calculations  based  on  ambient 
temperature,  solar  radiation,  wind  velocity,  ammunition  temperature,  magazine  and  pallet 
dimensions,  and  magazine  physical  property  data  (e.g.,  density,  specific  heat)  were  used  to 
determine  insulation  factors  between  the  ambient  environment  and  the  ammunition.  Once  the 
insulation  factors  are  determined,  ammunition  temperatures  can  be  predicted  solely  on  ambient 
temperatures,  wind  velocity,  and  solar  radiation  exposure.  Initially,  sine  functions  were 
developed  to  define  ambient  temperatures.  These  functions  relate  ambient  temperatures  to  the 
time  of  year  allowing  engineers  to  more  effectively  predict  temperatures  the  ammunition  will 
experience.  Using  this  information,  shelf  life  forecasts  can  be  calculated  analytically. 

Scientific  Basis 

The  travel  of  thermal  energy  through  materials  is  best  understood  by  the  laws  of  heat 
transfer.  The  study  of  heat  transfer  differs  from  thermodynamic  studies,  in  that;  it  focuses  on  the 
rate  of  energy  transmission,  whereas  thermodynamics  is  the  basic  science  that  deals  with  en¬ 
ergy,  matter,  and  their  interactions. 
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The  following  heat  transfer  concepts  were  considered: 

•  Heat  transfer  is  based  on  temperature  gradients  within  two  or  more  bodies 

•  Heat  energy  flows  from  a  high  temperature  region  to  a  low  temperature  region 
until  equilibrium  is  achieved 

To  predict  temperatures  in  an  ammunition  storage  facility  using  the  approaches  discussed 
in  this  study,  the  following  data  must  be  available: 

•  Ambient  temperature,  wind  velocity,  and  solar  radiation  data  for  a  given  loca¬ 
tion 

•  Magazine  dimensions  and  pallet  locations 

•  Physical  property  parameters  (e.g.,  specific  heat,  geometry,  and  density)  of 
magazine  makeup  and  its  contents 

The  temperature  of  ammunition  (even  if  stored  in  a  magazine)  is  affected  by  solar  radia¬ 
tion.  The  energy  coming  from  the  sun,  solar  energy,  reaches  the  earth  in  the  form  of  electro¬ 
magnetic  waves.  Upon  impact  with  a  material’s  surface,  solar  radiation  will  raise  the 
temperature  of  the  material  above  the  ambient  temperature  to  a  higher  radiation-induced  tem¬ 
perature. 


LOCATION  CHARACTERIZATION 

The  first  major  step  towards  classifying  a  magazine  storage  system  was  to  quantify  the 
environment  at  a  specific  magazine  location.  The  quantification  was  accomplished  by  develop¬ 
ing  models  that  express  the  ambient  temperature  as  a  function  of  the  Julian  time.  Tempera¬ 
tures  on  an  hourly  basis  were  desired  in  order  to  develop  a  true  equilibrium  relationship 
between  environmental  and  ammunition  storage  temperatures  for  a  given  geographic  location. 
To  account  for  seasonal  and  daily  temperature  changes  in  the  environment,  trigonometric  rnod- 
els  were  developed.  These  models  lend  themselves  to  better  integration  with  predictive  models 
and  allow  consolidation  of  data  into  summary  values  for  ease  of  computer  storage. 

Even  though  both  metric  and  English  units  are  used  throughout  this  study,  all  calculations 
were  accomplished  in  metric  units. 

Daily  Temperature  Cycle 

The  daily  temperature  cycle  model,  equation  1 ,  was  developed  as  a  general  equation  for 
any  location.  The  values  that  are  entered  into  the  model  are  used  to  develop  the  precise  equa¬ 
tion  for  a  specific  location.  The  cycle  begins  at  Julian  hour  1  (1  A.M.)  and  ends  at  Julian  hour  24 
(12  A.M.,  midnight).  However,  since  it  is  a  trigonometric  (sine  wave)  model,  the  cycle  repeats 
itself.  Therefore,  the  temperature  modeled  for  the  25th  hour  is  equivalent  to  that  of  the  1st  hour. 

’The  Julian  calendar  is  used  as  the  standard  for  overall  consistency.  Day  1  is  equivalent  to 
January  1st,  Day  32  is  February  1st,  and  leap  year  is  neglected  in  this  study.  Julian  hour  refers 
to  military  time,  where  Julian  hour  1  is  equivalent  to  0100  hrs. 
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T(HOUR ) = TAVEm,  +  Jz^iaas^yeiUsi 

^  |_  "hour 

where 

T(HOUR)  is  the  environmental  temperature  as  a  function  of  Julian  hour 
TAve  is  the  average  temperature  for  a  given  location 
RG  is  the  maximum  minus  minimum  temperature  range  for  a  given  location 

x  is  the  Julian  hour  (1  to  24) 

ALIGN  is  the  horizontal  centering  adjustment  to  Julian  calendar  base  on  geographic 
location  (note) 

n  is  the  number  of  data  points  (24  for  hourly  data,  365  for  daily  data) 

hour  all  variables  with  this  subscript  are  “as  a  function  of  hours" 

Note:  The  value  for  ALIGN  is  determined  by  minimizing  the  absolute  ERROR  between  the  ac¬ 
tual  temperature  and  the  model-derived  temperature  where 

ERROR  =  max  \Tact  -  TM0DEL  | 

Since  Tmodel  is  a  function  of  the  value  for  ALIGN,  each  ALIGN  value  will  correspond  to  one 
ERROR  value.  Therefore,  we  vary  the  value  lor  ALIGN  until  we  achieve  the  minimum  ERROR 
value.  ( ALIGNhour  will  be  between  0  and  24,  and  ALIGNDay  will  be  between  0  and  365.)  The 
ALIGN  value  with  the  minimum  ERROR  value  is  inserted  into  the  Tmodel  equation  (as  a  con¬ 
stant)  to  characterize  a  particular  data  set. 

As  an  example,  the  daily  temperature  cycle  for  Kuwait  City  was  modeled  using  equation  1 . 
Published  historical  data  for  the  hourly  temperatures  for  Kuwait  City  were  the  basis  of  example 
equation  1A,  which  is  the  precise  equation  for  this  example 

T(HOUR)  =  98.4 +  y sin  MX|-*  + 15)  (1A) 

Figure  1  depicts  the  accuracy  of  the  model,  by  plotting  the  actual  hourly  temperatures  and 
the  modeled  daily  temperature  cycle  for  Kuwait  City.  Although  the  model  does  not  exactly  corre¬ 
late  to  the  actual  data  points,  the  T(HOUR)  graphed  line  and  actual  temperature  data  are  al¬ 
ways  within  3°F.  It  should  be  noted  that  the  highest  temperatures  for  the  Kuwait  City  cycle  are 
during  midday  between  1400  and  1600  hrs  (2  to  4  P.M.). 
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Figure  1 

Hourly  temperature  cycle  for  Kuwait  City 


Yearly  Temperature  Cycle 

The  next  step  is  to  develop  a  model  for  the  yearly  or  (day  to  day)  temperature  cycle  for  a 
specific  geographic  location.  To  maintain  consistency,  this  model  is  also  a  function  of  the  Julian 
date.  Equation  2  is  a  modification  of  the  T(HOUR)  model,  with  hourly  temperatures  replaced  by 
daily  temperatures.  In  this  case,  the  temperature  for  each  day  is  modeled  and  the  cycle  consists 
of  365  days  as  opposed  to  24  hrs.  Day  1  is  January  1st  and  day  365  is  December  31st.  This 
model  also  repeats  itself,  but  instead  of  repeating  from  day  to  day,  it  repeats  from  year  to  year. 
Therefore,  day  366  has  the  equivalent  temperature  as  day  1 . 

T(DAY)  =  Tavs +  ”a«rs  in 

where 


2#  (Yday  AUGNMY ) 


n, 


DAY 


(2) 


T(DAY)  is  the  environmental  temperature  as  a  function  of  Julian  day 
y  is  the  Julian  day  (1  to  365) 

DAY  all  variables  with  this  subscript  are  “as  a  function  of  days” 

As  an  example,  the  same  location  (Kuwait  City)  was  used  to  create  a  precise  yearly  tem¬ 
perature  cycle.  In  this  model,  it  was  apparent  that  the  highest  temperatures  occur  in  the  summer 
months,  whereas,  for  the  daily  cycle  the  highest  temperatures  occurred  during  midday  (as  ex¬ 
pected).  Published  historical  data  for  the  temperature  in  Kuwait  City  was  also  used  to  generate 
example  equation  2A  for  the  yearly  cycle. 


T(DAY)  =  9  8  +  ysin 


2^~(yi-365  228) 


365 


(2A) 
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Figure  2  shows  the  yearly  environmental  temperature  cycle  for  Kuwait  City.  Actual  data  is 
not  depicted  for  clarity.  The  hottest  days  are  between  Julian  date  200  (June  19th)  and  250  (Sep¬ 
tember  7th). 


Daily  temperature  cycle  for  Kuwait  City 

Composite  Yearly  Temperature  Cycle 

A  temperature  model  now  exists  for  every  hour  within  a  day  (equation  1)  and  for  every  day 
within  a  year  (equation  2).  By  combining  equations  1  and  2,  a  location  can  be  temperature 
‘characterized’  on  an  hourly  basis  for  an  entire  year.  The  combined  daily  and  yearly  temperature 
cycles  result  in  equation  3,  in  which  the  T(DAY)  model  replaces  the  average  hourly  temperature 

(Tavehour)  in  the  T(HOUR)  model. 


T (DAY, HOUR)  =  T(DAY)+  ™UR  sin 


^■n  {XHOUR  +  ALIGNhqur  ) 
^HOUR 


(3) 


where 


T(DAY,HOUR)  is  the  environmental  temperature  as  a  function  of  Julian  day  and  hour. 

Continuing  with  the  Kuwait  City  example,  an  overall  yearly  temperature  cycle  was  created. 
In  this  example,  the  diurnal  temperature  variation  (hourly  during  a  day)  is  integrated  with  the  an¬ 
nual  temperature  variation  (daily  during  a  year)  to  provide  a  complete  cycle  that  accounts  for 
each  and  every  hour  for  a  year.  The  combined  cycles  create  one  large  sine  wave  from  Julian 
Day  1  to  Day  365.  Inside  this  large  sine  wave,  there  are  364  small  sine  waves  between  each 
and  every  Julian  day.  The  overall  model  runs  for  the  entire  year  and  repeats  from  year  to  year. 
The  following  precise  example,  equation  3A,  for  Kuwait  City  results  in  8,760  (24x365)  data 
points,  accounting  for  every  hour  within  a  year. 
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T (DAY, HOUR)  =  98.0  +  ^sin 


2tc  (yi_365  +  228) 

24 

+ — ! 

’2^(x,_24  +  15) 

365 

2 

24  J 

(3A) 


Figure  3  depicts  an  interval  (from  Julian  date  180  to  230)  of  the  example  equation  3A  plot 
for  the  combined  model. 
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Figure  3 

Daily  temperature  cycle  for  Kuwait  City 


MAGAZINE  CHARACTERIZATION 

The  second  major  step  towards  classifying  a  magazine  storage  system  was  to  quantify 
the  magazine  itself,  which  is  accomplished  through  model  development  of  the  physical  and  di¬ 
mensional  properties  of  the  magazine.  These  properties  are  used  to  determine  the  insulation 
factors  for  a  magazine,  which  are  in  turn  used  to  determine  the  temperatures  at  key  locations 
within  the  maga-zine.  The  remainder  of  this  study  will  be  based  on  the  igloo  magazine  dis¬ 
cussed  in  this  section,  which  is  located  in  Hawthorne,  Nevada. 

Physical  and  Dimensional  Properties  Example 

An  igloo  magazine  was  used  as  an  example  for  ‘Magazine  Characterization.’  Dimensional 
properties  were  obtained  from  the  Industrial  Operations  Command  (IOC),  1994,  Storage  Facili¬ 
ties  report.  The  physical  properties  were  obtained  from  chemical  handbooks  for  the  materials 
that  make  up  the  igloo  magazine  wall,  the  earth  that  covers  the  igloo,  the  air  between  the  wall 
and  the  ammunition  pallets,  and  the  ammunition  pallets.  Physical  and  dimensional  property  val¬ 
ues  are  summarized  in  appendix  A. 
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Figure  42  is  a  photograph  of  the  Hawthorne  earth-covered  igloo  magazine  that  will  be 
used  as  an  example  for  the  remainder  of  this  study.  This  figure  also  shows  the  exterior  weather 
station,  which  measured  ambient  temperature,  wind  velocity,  and  solar  radiation. 

The  igloo  magazine  dimensions  are  depicted  in  figure  5.  The  specifying  dimensions  are 
60  ft  8  in.  L  x  26  ft  6  in.  W  x  12  ft  9  in.  H  igloo  magazine,  where  H  refers  to  the  highest  point  of 
an  arched  roof. 

Figure  62  is  a  photograph  of  the  interior  of  the  earth-covered  igloo  magazine,  which  shows 
the  ammunition  pallet  stacking  and  the  temperature  instrumentation. 

Figure  72  shows  a  computerized  isometric  view  of  the  pallet  stacking  procedure. 

Figure  82  depicts  the  exact  temperature  instrumentation  locations  within  the  igloo  maga¬ 
zine. 


Figure  9  is  a  Pro/ENGINEER  schematic,  which  shows  the  heat  transfer  layers  within  the 
igloo  magazine.  The  bar  on  the  left  is  a  blow-up  of  a  very  narrow  center  column  of  the  igloo. 
Models  were  generated  to  determine  a  temperature  at  the  top  of  each  layer  for  the  center  col¬ 
umn. 


Figure  4 

Weather  station  positioned  on  top  of  igloo 


2Figures  4  and  6  through  8  were  obtained  from  the  Defense  Ammunition  Center  website 
(www.dac.army.mil). 
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Figure  9 

Igloo  heat  transfer  layers 


AMBIENT  CONDITIONS 


The  ambient  conditions  experienced  by  an  ammunition  storage  magazine  include  the  ef¬ 
fects  of  solar  radiation  and  wind  velocity  in  addition  to  the  ambient  temperature.  Thankfully,  this 
type  of  data  was  obtained  for  many  ammunition  storage  locations  throughout  the  world  every 
half-hour  for  entire  years.  This  data  can  be  very  useful  for  engineers  to  assess  the  storage  con¬ 
ditions  that  the  ammunition  is  facing  and  its  apparent  shelf  life.  However,  the  number  of  data 
measurements  for  1  yr  alone  is  over  17,000  for  each  effect.  This  number  of  data  points  can  be 
very  unwieldy  for  computer  storage  even  for  today’s  technology.  The  models  below  were  devel¬ 
oped  to  eliminate  the  need  for  large  databases  by  consolidating  the  data  into  equation  form.  The 
models  for  ambient  temperature,  solar  radiation,  and  wind  velocity  are  intended  to  be  general 
for  all  locations;  although,  precise  models  for  Hawthorne  are  also  provided  in  the  subsequent 
sections.  Actual  (Jan  1, 1996)  versus  model-generated  (Day  1)  temperatures,  solar  radiations, 
and  wind  velocities  are  given  in  appendix  B. 


Ambient  Temperature  Model 

The  yearly  composite  equation  (equation  4)  and  precise  model  (example  equation  4A)  for 
the  ambient  temperature  ( Tamb )  experienced  by  the  Hawthorne  igloo  magazine  are 


'  AMB 


=  Tt 


AVE 


rgday  . 

+ - 

2n{y  +  ALIGN DAY) 

+  RGhour  sin 

2u  (x  +  ALIGN HQUR ) 

2 

nDAY 

2 

nHOUR 

(4) 
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where 


Tamb 

Tave 

RG 

ALIGN 

n 

HOUR 

DAY 

X 

y 


is  the  ambient  temperature  as  experienced  in  Hawthorne 

is  the  average  daily  temperature  in  Hawthorne 

is  the  maximum  minus  minimum  temperature  range  for  Hawthorne 

is  the  horizontal  centering  adjustment  to  Julian  calendar  for  Hawthorne 

is  the  number  of  data  points  (24  for  hourly  data,  365  for  daily  data) 

all  variables  with  this  subscript  are  “as  a  function  of  hours” 

all  variables  with  this  subscript  are  “as  a  function  of  days” 

is  the  Julian  hour  (1  to  24) 

is  the  Julian  day  (1  to  365) 

2^(x+14) 


TAMB'=  55. 7923 +  17. 11 75  sin 


2?r(y +  270) 
58.0916 


+  11. 4450 sin 


3.8197 


Figure  10  graphically  depicts  example  equation  4A  for  an  entire  year. 


Figure  10 

Ambient  temperature  versus  Julian  date 


(4A) 
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Solar  Radiation  Model 


Heat  transfer  from  the  ambient  environment  to  the  earth  that  covers  the  igloo  magazine  is 
exposed  to  solar  radiation  and  wind  convection.  Solar  radiation  values  were  recorded,  and  a 
solar  radiation  sine  wave  model  (equation  5)  was  established  as 


SRamb  =  I^^day  | + SR  day  +  RGqay  (|SRhour  |  +  SRhour  +  L  OCA  TE ) 


(5) 


where 


SRamb  is  the  solar  radiation  model  experienced  by  Hawthorne 


SR  day  ~ 


2^(y+/tZ./GA/My) 


4/1 


DAY 


,  is  the  solar  radiation  as  a  function  of  days 


SRhour  ~  s'n 


2ff,(x + ALIGN HqUR  ) 

^  n HOUR 


,  is  the  solar  radiation  as  a  function  of  hours 


x  is  the  Julian  hour  (1  to  24) 

y  is  the  Julian  day  (1  to  365) 

ALIGN  is  the  horizontal  centering  adjustment  to  Julian  calendar  for  Hawthorne 

LOCATE  is  the  vertical  centering  adjustment  to  Julian  calendar 

The  precise  model  (example  equation  5A)  for  the  solar  radiation  model  {SRamb)  experi¬ 
enced  by  the  Hawthorne  igloo  magazine  is 


SRamb  ~ 


4  *+.?"> ) +4JttM.j+«Ls44 m +si»(^ 

v  232.3664  y  V232.3664J  [_  U.8197;  U. 


x+18^ 
.8197  J 


-2.1 


(5A) 


Note:  The  solar  radiation  (SRamb)  sine  wave  was  modified  to  provide  positive  values  only.  The 
following  algorithm  was  entered  into  the  equation  model  in  order  to  negate  nonexistent  negative 
solar  radiation  values. 


IF(SRamb  <  0 ,SRamb  =  0);ELSE(SRamb  =  SRAMB ) 


Figure  1 1  shows  both  the  actual  and  modeled  solar  radiation  data. 

Figure  12  is  a  blown-up  section  of  figure  11,  which  is  intended  to  show  the  detailed  daily 
curves. 
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Julian  Date  (Days) 
Figure  1 1 

Solar  radiation  versus  Julian  date 
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Figure  12 

Solar  radiation  versus  Julian  date  (section  of  fig.  11) 
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Wind  Convection  Model 

Wind  velocity  values  were  also  recorded,  and  a  wind  velocity  ( WVamb )  sine  wave  model 
(eq  6)  was  established  as 

WVMB  =RG+ ^a^sinf  2*(*mu*+AUGN*»*)'}  (6) 

2  L  nHOUR 


where 

WVamb  is  the  wind  velocity  model  experienced  by  Hawthorne 
RG  is  the  maximum  minus  minimum  wind  velocity  range  for  Hawthorne 

ALIGN  is  the  horizontal  centering  adjustment  to  Julian  calendar  for  Hawthorne 

hour  all  variables  with  this  subscript  are  “as  a  function  of  hours” 

DAY  all  variables  with  this  subscript  are  “as  a  function  of  days" 

n  is  the  number  of  data  points  (24  for  hourly  data) 


x  is  the  Julian  hour  (1  to  24) 

The  precise  wind  velocity  (WVamb)  model  (example  equation  6A)  experienced  by  the  Haw¬ 
thorne  igloo  magazine  is 


WVAUB  =  8.8060  +  8.1825sin(  4^ 
amb  V3.8197 


Figure  13  shows  both  the  actual  and  modeled  wind  velocity  data. 
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Julian  Date  (Days) 


Figure  13 

Wind  speed  (mph)  versus  Julian  date 

The  wind  velocity  model  was  used  to  determine  the  coefficient  of  heat  transfer  at  ambient 
conditions  ( hAMB ).  Table  1  provides  the  sample  calculations  for  air  properties  at  both  300  K  and 
250  K.  Even  though  hA mb  can  be  considered  a  function  of  temperature,  for  air  it  only  varies 
slightly  as  shown  by  the  following  example.  As  a  result,  subsequent  calculations  will  use  the 
physical  property  constants  of  air  at  300  K. 


Table  1 

Coefficient  of  heat  transfer  (convection)  as  a  function  of  wind  velocity 


Physical  properties  of  air  at  300  K  (80°F) 

Physical  properties  of  air  at  250  K  (-1 0°F) 

y  =  1.59  x  1 0' W/s 

v=  1.14  x  1 0'W/s 

k  =  2.63  x  1 0'2W/(mK) 

k  =  2.23  x  10'2W/(mK) 

Pr=  0.707 

Pr=  0.72 

where 


v  =  kinematic  velocity 
k  =  thermal  conductivity 

_ Pr~  Prandtl  number 

Wind  velocity 

v  =  1.8  m/s  (4  mph) _ , 

Surface  depth 

L  =  0.6096  m  (2  ft) _ 

Reynolds  number 

(must  be  less  than  5  x  105  v-L/v) 

Re  =  6.91  x  104 _ 

Nuselt  number 


Wind  velocity 

1.8  m/s  (4  mph) _ 

Surface  depth 

L  =  0.6096  m  (2  ft) _ 

Reynolds  number 

(must  be  less  than  5  x  105  v-L/v) 

Re  =  6.59  x  104 _ 

Nuselt  number 


(Nu  =  0.332  Re0  5-P/°333) 


(Nu  =  0.332 •Re°-S-P/’-333) 


Nu  =  77.7 _ 

Convection  coefficient 


Nu  =  77.7 _ 

Convection  coefficient 


^amb  -  Nu-kA.) 

hAMS  =  3.36  W/(m2-K) 


(J hAMB  =  Nuk/L) 

hwe  =  3.36  W/(m 


iKI 


HEAT  TRANSFER 


The  following  will  use  the  sine  models  previously  developed  for  the  ambient  conditions  to 
calculate  internal  magazine  temperatures.  The  heat  transfer  theory  outlines  the  procedure  by 
which  the  integrated  relationships  are  obtained.  The  derivation  is  based  on  the  theories  dis¬ 
cussed  in  the  publication,  Heat  Transmission  by  William  H.  McAdams  (McGraw  Hill,  1964).  The 
heat  transfer  approach  is  employed  for  all  material  layers  with  a  modification  used  for  the  heat 
transfer  between  ambient  conditions  and  the  earth  layer  of  the  igloo  magazine  to  account  for 
solar  radiation  and  wind  velocity  effects. 

Heat  Transfer  Theory 

t 

A  simple  case  of  unsteady  heat  transfer  is  discussed.  Consider  a  material  layer  of  volume 
V,  surface  area  A,  and  thickness  L  at  unknown  temperature  Tin  contact  with  another  material 
layer  at  temperature  T2.  At  any  time  f,  the  quantity  of  heat,  dQ,  transferred  in  the  short  time,  dt, 
depends  upon  the  surface  area  of  the  first  material  layer,  the  difference  in  temperature  between 
the  two  layers  (T2  -  T),  and  the  coefficient  of  heat  transfer  between  the  two  layers  (h).  There¬ 
fore,  by  determining  a  heat  balance  on  the  first  layer,  with  a  density/? and  specific  heat  c yields: 


Step  1:  dQ  =  h-A(T2-T)dt  =  VpcdT 

Assuming  h,  A,  V,  p,  and  c  are  constant,  define  U  (the  overall  heat  transfer  coefficient): 
hA 

Step  2:  U  =  - — 

Vpc 

Substitute  step  2  into  step  1  and  rearrange: 

Step  3:  Udt  =  -r-  1  -rdT 


Integrate  from  t  =  0  to  t  =  V  and  T=  T  to  7=  7/ 

s,ep4:  KfeV 

Integration  yields: 

Step  5:  Ut'  =  lnLzL' 

'2  7, 

Solving  for  7/  gives: 

Step  6:  77=  7", 

Where  7/  is  the  surface  temperature  of  the  first  material  layer  at  time  t' . 
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Material  Layer  Heat  Balances 

The  heat  balance  (equations  7  to  1 1)  is  the  starting  equation  in  the  determination  of  inte¬ 
rior  magazine  temperatures.  Balances  were  performed  around  the  outer  surface  of  the  earth, 
wall,  interior  air,  outer  pallet,  and  inner  pallet  layers.  The  equation  for  heat  transfer  between 
ambient  conditions  and  the  earth  layer  has  an  additional  expression  (cceAambSRamb)  to  account 
for  solar  radiation  effects.  All  other  material  layer  derivations  are  based  on  the  heat  transfer  the¬ 


ory. 

=  \_^AMB^AMB  ( TaMB  ~  J  dt  =  ^AMbPaMB^AMB^  (^) 

dQ^„=heAE(TE-Tw)dt=VEpEcedT  (8) 

dQw_+A  =  fyyAv  0~w  "■  TA)dt  =  VwpwcwdT  (9) 

dQA^0  =  hAAA  (Ta  -T0)dt  =  VApAcAdT  (10) 

dQ0->i  =  h0Ac  (T0  -T,)dt=  V0p0c0dT  (11) 


where 

x  is  the  subscript  referring  to  material  layer: 

x  =  AMB  for  ambient  conditions 
x  =  E  for  earth 
x=W  for  wall 

x  =  A  for  air  (inside  the  igloo  magazine) 
x  =  O  for  the  outer  pallet  layer 
x  =  /  for  the  inner  pallet  layer 

hx  is  the  coefficient  of  heat  transfer  of  material  layer  x  [W/(m2  K)] 

Ax  is  the  area  of  material  layer  x  through  which  heat  flows  at  right  angles  (m2) 

ax  is  the  solar  absorptivity  of  the  earth  layer  (unitless) 

Vx  is  the  volume  of  material  layer  x  (m3) 
is  the  density  of  material  layer  x  (kg/m3) 
cx  is  the  heat  capacity  of  material  layer  x  [J/(kg-K)] 

Overall  Heat  Transfer  Coefficients 

Once  the  properties  and  dimensions  are  ascertained,  the  overall  heat  transfer  coefficients 
(U)  can  be  calculated  (eqs  12  to  16)  for  each  heat  transfer  layer.  Note  that  for  solid  layers  the 
coefficient  of  heat  transfer  (h)  is  equal  to  the  material  heat  conductivity  {k)  divided  by  the  thick¬ 
ness  of  the  layer  (L),  h  =  k/L. 
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(12) 


1 I  _  ^AMB^AMB 
UAMB  ~  • 


VambPambcamb 


kgA 

LeVbPbc 


(13) 


L-W^W  Pw^W 

(14) 

u 

*  LaVaPaCa 

(15) 

y  k0Ao 

LcYoPoco 

(16) 

where 

Ux  is  the  overall  heat  transfer  coefficient  of  material  layer  x  (sec-1)  (table  2) 

kx  is  the  thermal  conductivity  of  material  layer  x  [W/(mK)]  . 

Lx  is  the  thickness  of  material  layer  x  (m) 

Physical  and  dimensional  property  values  are  in  appendix  A.  By  substituting  values  into 
equations  12  through  16,  the  overall  heat  transfer  coefficients  are  determined  for  the  igloo 
magazine. 

Table  2 

Overall  heat  transfer  coefficients  for  the  material  layers 

Uaub  1  IjhZS.  =  1  -85x10-*  to  9.67x1 0-a  sec-1 
l/E=  3.65x1 0~7  sec-1 
1.87x10*  sec-1 
1/4=  1.1 2x1  (p- 'sec-1 
Lfo=  1.67x10^  sec-1 

Igloo  Magazine  Temperature  at  Different  Locations 

The  overall  heat  transfer  coefficients  ( U)  allow  the  calculation  of  the  heat  transfer  through 
each  material  layer,  enabling  the  calculation  of  temperatures  at  various  locations  inside  the 
magazine  as  a  function  of  the  ambient  temperature,  solar  radiation,  and  wind  speed.  Values  for 
Julian  day  1  are  listed  in  appendix  B.  The  temperatures  were  calculated  using  equations  17 
through  21  which  represent  the  culmination  of  the  preceding  heat  transfer  derivations.  Model¬ 
generated  temperatures  at  half-hour  intervals  for  Julian  day  1  for  the  five  locations  are  provided 
in  appendix  C. 
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By  integrating  equations  7  through  1 1  depicted  in  Material  Layer  Heat  Balances,  equa¬ 
tions  17  through  21  are  derived.  The  derivation  of  equation  17  is  shown  in  appendix  D. 


where 

x  is  the  subscript  referring  to  material  layer: 

x  =  AMB  for  ambient  conditions 
x  =  E  for  earth 
x=  W for  wall 

x  =  A  for  air  (inside  the  igloo  magazine) 
x  =  O  for  the  outer  pallet  layer 
x  =  /  for  the  inner  pallet  layer 

Tx  is  the  temperature  of  material  layer  x  at  time  t'  (°F) 

Tx  is  the  temperature  of  material  layer  x  at  time  =  0  (°F) 

Ux  is  the  overall  heat  transfer  coefficient  of  material  layer  x  (sec-1) 
t'  is  the  selected  time  of  temperature  calculation  (sec) 


T>E  =  TEe 


-UampF 


Tamb  + 


aB l(1_e-(W') 

r,AMB 


where 


(17) 


T'e  is  the  outer  surface  temperature  of  earth  layer  at  time  t'  (°F) 

Te  is  the  outer  surface  temperature  of  earth  layer  at  time  =  0  (°F) 

Tamb  is  the  ambient  temperature  at  time  =  0  (°F) 

aE  is  the  solar  absorptivity  of  the  earth  layer  =  0.1 0  (unitless) 

SRamb  is  the  observed  solar  radiation  (W/m2) 


T'  -T 

T’  -T  E  'w 
1 w-'e  euev 


(18) 


where 

is  the  outer  surface  temperature  of  wall  layer  at  time  t'  (°F) 
Tw  is  the  outer  surface  temperature  of  wall  layer  at  time  =  0  (°F) 
T  -T 

T<  _  T  w  A 

1 A  ~  'W  eUwf 


(19) 
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where 


Ta  is  the  outer  surface  temperature  of  air  layer  at  time  t'  (°F) 

Ta  is  the  outer  surface  temperature  of  air  layer  at  time  =  0  (°F) 

T'  -T 

T>  -T  A  0 

'o  -  A  eUAV 


(20) 


where 

To  is  the  outer  surface  temperature  of  outer  pallet  layer  at  time  t'  (°F) 
T0  is  the  outer  surface  temperature  of  outer  pallet  layer  at  time  =  0  (°F) 


T'  -T 
T'-T  0  1 

'/-,o  euQr 


(21) 


where 

7/  is  the  outer  surface  temperature  of  inner  pallet  layer  at  time  t'  (°F) 

7/  is  the  outer  surface  temperature  of  inner  pallet  layer  at  time  =  0  (°F) 

The  measured  temperature  data  at  external,  earth-covered,  and  internal  (middle  pallet) 
locations  for  the  Hawthorne  igloo  magazine  for  a  year  are  graphically  displayed  in  figure  14. 
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Figure  14 

Actual  temperature  versus  Julian  date 
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The  application  of  equations  17  through  21  to  the  Hawthorne  igloo  example  are  depicted 
in  figure  15,  as  modification  of  the  previous  Pro/ENGINEER  schematic  (fig.  9).  Calculated  val¬ 
ues  for  ambient  conditions  and  for  each  layer  were  included  for  a  randomly  selected  point  in 
time  (Julian  day  305,  hour  21). 


Figure  15 

Igloo  heat  transfer  layers  (Julian  day  305,  hour  21) 


Figure  16  graphically  compares  the  7/  model-generated  temperatures  and  the  middle 

pallet  location  temperature  measurements.  Figure  16  shows  only  the  end  results  of  the  transfer 
model  (eq  21),  which  is  dependent  upon  all  previous  equations.  Note  that  the  results  obtained 
for  the  7/  model  are  independent  of  the  results  from  the  middle  pallet  measurements.  The  77 
model  was  developed  based  on  the  transfer  models  and  the  physical  property  data  associated 
with  the  igloo.  Whereas,  the  middle  pallet  location  values  are  thermo-couple  measurements  that 
were  obtained  from  the  DACS  website  (www.dac.army.mil). 


Figure  16 

Pallet  temperature  (actual  versus  model)  versus  Julian  date 
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Figure  17  is  a  histogram  that  compares  the  7/  model  temperature  output  to  the  measured 
middle  pallet  layer  temperature  from  figure  14.  The  average  difference  between  the  two  curves 
is  2.08°F  and  the  maximum  difference  is  5.60°F.  Note  that  the  values  in  the  histogram  were 
rounded  to  the  nearest  whole  number. 

30%  1— — — . . . . 1 


Temperature  Difference  (F) 
Figure  17 

Temperature  difference  histogram 


PREDICTIVE  ENGINEERING  METHODS 

The  Arrhenius  equation  is  used  to  determine  an  item’s  degradation  rate.  The  Life  Con¬ 
sumed  Concept  calculates  the  amount  of  shelf  life  consumed  for  a  given  storage  environment, 
using  information  from  the  Arrhenius  equation.  Two  examples  are  provided  to  demonstrate  the 
use  of  the  Arrhenius  equation  and  the  Life  Consumed  Concept. 

Arrhenius  Equation 

The  Arrhenius  equation  (eq  22)  is  used  to  model  the  dependence  of  an  item’s  degradation 
rate  constant  ( k )  on  temperature  (7).  This  equation  is  graphed  as  a  straight  line  when  the  loga¬ 
rithm  of  the  rate  constant  is  plotted  against  the  reciprocal  of  the  absolute  temperature.  The  pre¬ 
exponential  factor  (A)  and  the  activation  energy  ( Ea )  are  calculated  based  on  experimentally  de¬ 
termined  rate  constant/temperature  sets.  Once  A  and  Ea  are  determined  for  a  given  item,  the 
degradation  rate  constant  can  be  determined  for  any  temperature.  Once  the  rate  constants  are 
known,  the  life  prediction  can  be  calculated. 

k  =  A-e^  (22) 
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where 


k  is  the  degradation  rate  constant  for  a  given  item  (percent  of  life  consumed  per  year) 

A  is  the  pre-exponential  factor  (percent  of  life  consumed  per  year) 

Ea  is  the  activation  energy  [BTU/(lbmol)] 

R  is  the  universal  gas  constant  [1 .9858  BTU/(lb  mol  °R)] 

T  is  the  exposure  temperature  (°R) 

Temperatures  (7)  with  their  corresponding  degradation  rate  constants  (k)  are  inserted  into 
the  Arrhenius  equation  to  determine  A  and  Ea.  With  a  minimum  of  two  ( T  with  k)  data  sets,  the  A 
and  E„  are  calculated  using  two  equations  with  two  unknowns,  by  solving  equation  22  for  each 
data  set. 

Example:  Given  degradation  rates  (k)  for  an  ammunition  item: 

keo  =1/0.15  (percent  of  life  consumed/year)  at  80°F  (corresponds  to  a  15-yr  life) 
kgo  =  1/0.05  (percent  of  life  consumed/year)  at  906F  (corresponds  to  a  5-yr  life) 

Ea  =  64,716  BTU/(lbmol) 

A  =  1.1 2x1 027  (percent  of  life  consumed/year) 

Note:  The  keo  and  /c90  values  are  hypothetical  and  selected  for  this  conceptual  example  only. 
Actual  values  would  be  determined  through  an  accelerated  aging  experiment.  An  accelerated 
aging  experiment  is  performed  on  an  item  by  placing  it  in  a  high  stress  (temperature)  chamber 
for  a  predetermined  period  of  time.  This  process  is  designed  to  simulate  long-term  storage  ef¬ 
fects  using  high  stress  short-term  conditioning. 

Life  Consumed  Concept 

The  Ea  and  A  values  are  entered  into  equation  23  to  predict  the  overall  life  consumed.  The 
life  consumed  ( Lt )  is  the  percent  of  shelf  life  that  a  particular  item  has  used.  Note  that  equation 
23  assumes  that  the  item  follows  a  zero-order  degradation  rate. 


In 


=  kt  +  L0=Ae 


t +Lq 


(23) 


23 


where 


L,  is  the  total  percent  of  life  consumed  (percent) 
t  is  the  exposure  time  (years) 

L0  is  the  initial  or  previous  percent  of  life  consumed  (percent) 

Solving  equation  23  for  the  exposure  time  ( t )  provides  (equation  24)  remaining  life  for  a 
given  temperature  environment. 


L,  -Z-0  L,-4 

k  Ae^ 


(24) 


This  approach  allows  calculations  to  be  obtained  for  different  degradation  rates  and  expo¬ 
sure  temperatures. 

Example:  If  a  newly  produced  item  is  stored  for  2  yrs  in  a  90°F  environment  and 
then  stored  at  80°F,  determine  its  remaining  life  for  the  80°F  environment  by  using  equations  23 
and  24. 


Stepl  :  L(i9o  =  /Cgofgo +  Lo  (eq  23) 


1 

0.05 


•2  +  0 


Lf, go  =  40.0%  of  life  consumed  at  90°F 


Step  2.  f80  — 


^i,8o  h,< 


.90 


'80 


(eq  24) 


ho  ~ 


100-40.0 

1/ 

70.15 


Teo  =  9  yrs  of  life  remaining  at  80°F 


If  a  newly  produced  item  with  the  aforementioned  degradation  rate  constants  (k60  and  k90) 
experiences  a  2-yr  deployment  in  a  90°F  environment,  it  will  survive  nine  additional  years  when 
subjected  to  an  80°F  environment.  Figure  18  graphically  depicts  the  solution  for  the  previous 
example  and  provides  plots  for  the  cases  in  which  the  ammunition  is  solely  stored  in  90°F  or 
80°F  environment. 
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Exposure  Time  (year^) 

Figure  18 

Life  consumed  versus  exposure  time 


As  a  second  example,  the  Life  Consumed  Concept  can  be  applied  to  calculate  the  differ¬ 
ence  in  life  consumed  for  items  stored  at  Hawthorne  at  the  two  different  temperature  profiles. 
For  this  example,  the  two  temperature  profiles  are  ambient  temperature  (TAmb)  (i.e.,  outside  un¬ 
protected  not  in  direct  sunlight)  and  the  outer  pallet  inside  an  igloo  temperature  (Tq ).  Based  on 
the  hypothetical  Arrhenius  rates  previously  calculated,  the  percent  life  consumed  as  a  function 
of  exposure  time  for  the  two  temperature  profiles  is  graphically  presented  in  figure  19.  Note  that 
both  plots  are  curves  due  to  the  changing  temperature  over  the  course  of  a  year.  The  cumula¬ 
tive  annual  effect  at  TAMB  is  2.69%  life  consumed  and  at  Tq  is  2.56%  life  consumed.  By  using 
equation  23,  the  TAMB  exposed  item  will  have  a  37-yr  life  and  the  Tq  exposed  item  will  have  a 
39-yr  life.  For  this  example,  the  2-yr  difference  in  life  is  solely  due  to  the  different  temperature 
profiles.  The  effects  of  solar  radiation  (which  can  be  substantial),  precipitation,  and  other  envi¬ 
ronmental  conditions  (e.g.,  blowing  sand  and  dust)  would  compound  the  degradation  rate  of  an 
unprotected  item. 


Exposure  Time  (years) 


Figure  19 

Life  consumed  in  1  yr 
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SUMMARY 


By  mathematically  characterizing  the  (1 )  yearly  ambient  temperature,  solar  radiation  and 
wind  convection  and  (2)  heat  transfer  through  the  layers  of  a  storage  magazine  (igloo),  a 
method  was  created  to  calculate  a  temperature  profile  for  igloo  stored  materiel.  Using  the  Ar¬ 
rhenius  equation,  a  Life  Consumed  Concept  was  created  to  predict  shelf  life  of  an  item  as  a 
function  of  its  temperature.  By  combining  these  two  analytical  methods,  an  item’s  shelf  life  can 
be  calculated  based  on  ambient  conditions  and  its  storage  facility. 

The  procedural  notions  presented  in  this  study  can  be  extended  to  include  an  economic 
analysis.  An  economic  analysis  would  focus  on  optimization  by  basing  logistic  decisions  on 
knowing  an  items  shelf  life. 
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APPENDIX  A 

PHYSICAL  AND  DIMENSIONAL  PROPERTY  VALUES 
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Parameter 

Ambient 

Condition 

Earth 

Wall 

Air 

Outer  Pallet 
Layer 

Le  =  0.61 

Lw=  0.61 

Z-o  =  0.76 

t&pflpgfl 

Aamb  =  0.37 

Ae  =  0.37 

Aw  -  0.37 

Aa  =  0.37 

A0  =  0.37 

Vamb  =  0.23 

VE  =  0.37 

VW  =  0.23 

V,,  =  0.28 

Vo  =  0.28 

Conductivity 

[W/(mK)] 

kE  =  0.27 

/%  =  1.40 

/c0  =  237.00 

Heat  capacity 

[J/kg-K)] 

Camb  =  1007 

Cjp  =  800 

c^v  =  880 

cA  =  1007 

c0  =  903 

Density  (kg/m3) 

Pamb  -  1-16 

R25ESEEIH 

ilWLkkl»l»TiW 

Pa  =  1.16 

Po  =  2702.0 
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APPENDIX  B 

ACTUAL  (1  JAN  1996)  VERSUS  MODEL-GENERATED  (DAY  1)  TEMPERATURE,  SOLAR 

RADIATIONS  AND  WIND  VELOCITIES 
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Julian  Date 


and  Time 


Scaled  Time 
(days) 


Temperature  (°F 


Sine  Wave 
Model 


Solar  Radiation 


Sine  Wave 


mi 


Ambient 


0.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


0.0000 


1.0000 


0.0000 


2.0000 


10.0000 


41.0000 


186.0000 


264.0000 


328.0000 


397.0000 


438.0000 


480.0000 


500.0000 


501.0000 


513.0000 


482.0000 


474.0000 


402.0000 


315.0000 


255.0000 


211.0000 


110.0000 


23.0000 


5.0000 


1.0000 


1.0000 


0.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


Model 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


103.4846 


213.5462 


310.0677 


391.3978 


456.1448 


503.2008 


531.7608 


541.3360 


531.7627 


503.2046 


456.1504 


391.4052 


310.0768 


213.5567 


103.4964 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


Wind  Speed  (mph 


Sine  Wave 
Ambient  Model 


4.26 


3.99 


4.33 


4.60 


3.61 


3.81 


5.38 


0.0000 


6.28 

15.89 

4.66 

16.37 

6.32 

16.71 

6.03 

16.92 

5.06 

16.99 

3.83 

16.92 

6.12 

16.71 

4.19 

16.37 

5.11 

15.89 

3.85 

15.30 

3.53 

14.59 

2.45 

13.79 

4.11 

12.90 

6.18 

11.94 

4.42 

10.92 

4.50 

9.87 

3.62 

8.81 

7.15 

7.74 

4.50 

6.69 

4.46 

5.67 

4.25 

4.71 

3.05 

3.82 
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APPENDIX  C 

MODEL-GENERATED  (DAY  1)  TEMPERATURES  FOR  HAWTHORNE  IGLOO  LAYERS 


37 


Julian  Date  and  Time 


Calculated  Surface  Temperatures 
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APPENDIX  D 

DERIVATION  OF  EQUATION  17  FROM  EQUATION  7 
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dQAMB-+E  “  \_^AMB^AMB  (j~AMB  +  aE  Aamb^^amb]^1 

Substitute  for  clarification: 

A  =  hAMB  •  Aamb 
Yamb  =  Tamb 

y- te 

B=  cce-  Aamb  *  SRamb 
dx=dt 

C  =  Vamb  •  Pamb  •  Camb 
Step  1 :  [A(yAUB  -y)+B]dx  =  Cdy 
Rearrange  and  set  equal  to  0 

Substitute: 

d=4 


_a  _e 
t_  Cy/WB  c 


Step  3:  —+Dy+E  =  0 
dx 

Rearrange  to  the  form  of  M  •  dx  +  N  •  ofy  =  0: 
Step  4:  (D  y  +  E)dx  +  dy  =  0 

In  this  form,  the  integration  factor  (//)  is  equal  to: 

M  =  Jf(x)dx 


Where: 


8M_+dN_ 
dy  dx 


-  ^ambPambcamb^  (?) 
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Example:  f(x)  =  y(D  +  0  )  =  D 

Therefore,  multiply  through  by  p  =  eDx : 

Step  5:  e°x (D-y  +  E)dx  +  eDxdy  =  0 

Rearrange  to  perform  the  Reverse  Chain  Rule  on  the  left-hand  side: 

Step  6:  eDx-+D-y-eDx  =  -E-eDx 
dx 

Perform  the  Reverse  Chain  Rule: 


Step 


7:  d(ye°‘).  £.e°- 

dx 

Rearrange  and  integrate  from  y  =  yo  to  y  =  yi  and  from  x  =  x0  to  x  =  Xi: 


S«eP8:  J^  =  -£K 


yo 


Note:  eDx°  =  1 ,  since  x0  =  0: 


Step  9:  y,  -eDx  -y0  =  (e0x  -l) 
Solve  for  yy. 


Step  10:  y1=-|(i_eDx)  +  yo  e0x 

Substitute  back  to  heat  transfer  form: 


^=y, 

Te  =  yb 
Uamb  -  D 
f'  =  x, 


7S 


r£e"uw'  + 


+ 


aE^^AMB 
^AMB  J 


(l-e_tw) 


(17) 
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